Weakly coupled one-dimensional Mott insulators 
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We consider a model of one-dimensional Mott insulators coupled by a weak interchain tunnelling 
t± . We first determine the single-particle Green's function of a single chain by exact field-theoretical 
methods and then take the tunnelling into account by means of a Random Phase Approximation 
(RPA). In order to embed this approximation into a well-defined expansion with a small parameter, 
the Fourier transform T±{k) of the interchain coupling is assumed to have a small support in mo- 
mentum space such that every integration over transverse wave vector yields a small factor acq ^ 1. 
When Tl (0) exceeds a critical value, a small Fermi surface develops in the form of electron and hole 
pockets. We demonstrate that Luttinger's theorem holds both in the insulating and in the metallic 
phases. We find that the quasi-particle residue Z increases very fast through the transition and 
quickly reaches a value of about 0.4 — 0.6. The metallic state close to the transition retains many 
features of the one- dimensional system in the form of strong incoherent continua. 

PACS numbers: 71.10.Pm, 72.80.Sk 

I. INTRODUCTION 



Ti(k) = -rx(k + Q) 



(1) 



The problem of the Mott metal-insulator transition, 
which is a consequence of electron interactions and not 
due to the bandstructure, has been the subject of great 
interestEl since the pioneering works by Mottlj. In 1964 
Hubbard suggested a solution of this problem for a par- 
ticular modelEl. The key point of Hubbard's approach was 
to take into account the on-site Coulomb interaction U 
in zeroeth order and develop perturbation theory in t/U. 
Unfortunately, the transition from the metallic to the in- 
sulating phase is expected to occur when the tunnelling 
matrix element t is comparable to the on-site repulsion 
U. In this regime expansions in either t/U ot U/t are 
not applicable. To overcome this difficulty, various ap- 
proaches have been suggested. One of them considers a 
SU(N)-symmetric generalization fttithe electronic model 
and carries out a l/A^-expansiorDiJ. Another approach 
is based on the so-called Dynamical Mean Field theory, 
which considers a lattice in infinitely many jdimensions 
D ^ oo (seeQ and references therein and alsoQH). 

In this paper we follow a different route. We consider 
a quasi-one-dimensional model of interacting electrons, 
where the tunneling along one direction is much larger 
than in all others. This model can be described in terms 
of weakly coupled chains. When the band is half-filled 
and the interchain tunneling is switched off, Umklapp 
processes dynamically generate a spectral gap M. The 
same mechanism can generate gaps at any commensurate 
filling e.g. quarter-filling, but only if the interactions are 
sufficiently strong. In the half-filled case the decoupled 
chains develop a Mott-Hubbard gap for any positive value 
oi U > 0. In what follows we shall assume that the 
interchain tunneling is weak but long-ranged, such that 
the Fourier transform of the interchain tunneling matrix 
element has two strong peaks in the Brillouin zone, one 
around zero wave vector and the other one around Q = 
(0, TT, tt) such that 



Due to the long-range character of the tunneling, the 
width of these peaks, denoted by kq, is small (kq ^ 
1). Hence every integration over the transverse mo- 
mentum yields a small parameter Kq. The main idea 
of our approach is to treat the individual chains non- 
perturbatively and then to employ the Random Phase 
Approximation (RPA) to take into account the interchain 
tunneling. Since the Mott metal-insulator transition is 
associated with development of a coherent single-particle 
excitation branch, the correlation function we are inter- 
ested in is the single-electron Green's function G. The 
RPA expression for G is 



G{u;,q,k) = [G^\LO,q)-T^{k)Y 



(2) 



where q is the momentum along the chain direction and 
Go is the single-particle Green's function for an individ- 
ual chain. The corrections to RPA are of higher order in 
Kg, which is the small parameter in our expansion. The 
smallness of Kq suppresses multiparticle tunneling pro- 
cesses, which generate exchange interactions between the 
chains and eventually lead to a three-dimensional phase 
transition. However, for small Kq the temperature at 
which this transition occurs is much smaller than the 
characteristic energy scale M of the problem. Therefore 
throughout the paper we shall assume that we work at 
temperatures much smaller than M (such that the ther- 
mal effects on the single-particle Green's function can be 
neglected) but much larger than the transition tempera- 
ture Tc, i.e. Tc <^T <^ M. We note that an analogous 
RPA has been used by Wencl to. study weakly coupled 
Luttinger liquids (LL), see alsoEj for a derivation based 
on functional integrals. An improved calculation for the 
case of coupled LLs has recently been carried out by 
ArrigoniO. It is based on taking into account multipoint 
correlation functions of upcoupled chains in a controlled 
way. In a Mott insulator multipoint correlation functions 
are much more difficult to determine than in a LL. We 
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therefore postpone the discussion of corrections to RPA 
until future pubhcations. 

We also neglect the long distance tails of the Coulomb 
interaction. This can be justified at finite temperatures 
if the dielectric constant at small frequencies is large. 

The basic input of our approach is Gq. To calculate 
this function we assume that for an isolated chain the 
Mott-Hubbard gap M is much smaller than the band 
width In this case one can use the continuous field 
theoretical description in the low energy limit. The corre- 
sponding field theory will described in detail in the next 
section and is highly universal; all information about the 
underlying lattice model is incorporated in just three di- 
mensionless constants. Thus, our assumptions can be 
summarized by the inequalities 

D\i > M{U) w ti(0) > t^{\k\ > Ko), 
T^^T<^ M (3) 

where ([/) symbolizes dependence of the Mott-Hubbard 
gap from the interactions. The smallness of M does not 
imply that the interactions are weak. For example, in 
the Hubbard model the field theoretical description wades 
well even for U/t\\ as large as 2-3 (see, for example,ll3). 
At such interaction strength one should expect the spin 
and charge velocities to be considerably different. This is 
a consequences of spin-charge separation, which is one of 
the most interesting features of one-dimensional strongly 
correlated systems. The effects of spin-charge separa- 
tion in the insulating regime deserve a special discussion, 
which will be given later in the paper. 

The outline of this paper is as follows. In section |l| we 
introduce the model and discuss the limit of decoupled 
chains. In section^ we take the interchain coupling into 
account and dete rmine various physics properties. 

In section |VIl| we discuss possible experimental appli- 
cations. 



II. THE MODEL 

We take as a starting point the following lattice model 
of correlated electrons 



(4) 



Here l,m label Hubbard chains and n labels the sites 
along a given chain. In the physically most interesting 
situation the interchain hopping matrix elements Um are 
taken to be equal to t±^ if I and m are nearest neighbours 
and zero otherwise. For this choice of interchain hopping 
our calculational scheme is uncontrolled: there is no small 
expansion parameter. One may nevertheless apply our 
scheme and hope that the corrections are small in some 
regime of temperatures. On the other hand one may 



choose the t„m's in such a way that a small expansion 
parameter is introduced, as explained above and in more 
detail in Appendix 

Let us first consider the case of uncoupled chains t_L — 
0. For weak repulsion U < t and low energies a field- 
theory description is appropriate. Keeping only modes 
in the vicinity of the Fermi momenta ±kp = ±7r/2ao, we 
may decompose the lattice Fermi operators as 



.(0 



exp 



(5) 



where ag is the lattice spacing. Inserting this prescription 
into the Hamiltonian and dropping the chain index 
{I) for the time being, one obtains 



-— I dx 
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J- J :] 



-1-25 J [I • I - J • J] 



(6) 



where vp = 2tao is the Fermi velocity and g = Uag. 
Here J and I are the chiral components of SU(2) spin 
and pseudospin currents 



-E 

2 ^ 

CT 



R\Rn 



/+ = r\r\ 



-R\Ri 



J+ -- 
J+ 



= R\Ri 



(7) 



By ernploying the Sugawara construction, the Hamilto- 
nian (H) can now be split into two parts, earresponding 
to the spin and charge sectors respectivelyli^ 



Hs = 



2ttv 



3 

27rUr 



dx \ : J ■ 3 : 



3 3 



— 2g J dx 3 ■ 3 , 
J dx[:I-I: + ■.i-i-]+2g J dxl-i. (8) 



Here — vp — C/ao/27r and Vc — vp + C/ao/27r. Apart 
from the (marginally) irrelevant current-current interac- 
tion in the spin sector and the difference in spin and 
charge velocities, the Hamiltonian (^ is identical to the 
one of the SU(2) Thirring model. The SU(2) Thirring 
model can be bosonized in terms of a Sine-Gordon model 
and a free boson; the resulting action density is 



^^^^ [v7\dr^s)^+Vs{d,^sf] , 
iuTT 



(9) 
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where [3 = 1. If we consider additional small density- 
density interactions between nearest neighbour sites in 
the field theory limit is again of the form but 
now (3 < 1. Using the integrability of the model ^ it 
is possible to determine dynamical correlation functions. 
This is the subject of the following section. 



III. UNCOUPLED CHAINS 

The calculation of the spectral function for half-filled 
Mott insulators is based on the following principles: 

(a) Locality: Ra,La are local fields. This allows us 
to employ the standard formfactor approachtfl. Some 
important elements of this approach are reviewed in Ap- 
pendix ^ 

(b) Spin-Charge separation: The Hamiltonian (||) is 
the sum of two parts representing the spin and charge de- 
grees of freedom respectively and by means of bosoniza- 
tion one can represent Ra,La as products of fields be- 
longing to the different sectors. For example: 



^^^'-^ e^P U<^c) exp 



(10) 



where 4'cT<t>s are chiral components of the bosonic fields, 
rjcr are Klein factors and the plus (minus) sign corre- 
sponds to (7 =] (a =i). The components into which 
the creation and annihilation operators are factorized are 
nonlocal fields, but the factorization of the field means 
that all its formfactors also factorize, at least in the limit 
when one of the sectors becomes massless. 

(c) "Triviality" of the spin sector: In the limit — > 
the field theory for the spin sector becomes massless and 
the correlation function of the spin exponent in (10) is 
simply given by 



exp 



exp 



1 



(11) 



(d) Lorentz invariance and Watson's theorem : Let us 
consider formfactors of the left moving fermi operator 
L„{x). The first nonvanishing formfactor is between the 
vacuum and a scattering state of one spinon (with spin a) 
and one antiholon. We denote the rapidity of the spinon 
(antiholon)by 6s {0c) (our notations are summarised in 
Appendix g). Lorentz invariance and Watson's theorem 
impose the following form for the first nonvanishing form- 
factor for the fcrmion operator 



exp[ 



.)/4] f{Oc 



(12) 



The function f(9) is periodic with period 2i7r and docs 
not coaJtain poles. Furthermore the usual asymptotic 
boundij yields 



lim (0\L\9r 



exp 



< const 



(13) 



where a = c, s. This leaves us with the only possibil- 
ity f{d) — const. As we see, this matrix element does 
not depend on the anisotropy of the interaction (up to 
a constant prefactor). The reason for this is that all in- 
formation about the anisotropy of the coupling constants 
is contained in the two-particle S-matrices of holons and 
spinous and does not influence the single-particle emis- 
sion. Aiuexplicit expression for / has recently been ob- 
tained inlij; for /J'^ — > 1 (the isotropic case) the numerical 
value is 



/ = (Zo/2^)i/2; Zo = 0.9218. 



(14) 



The matrix element (|12|) describes the emission of one 
kink in the charge and one kink in the spin sector. The 
charge part of this first term in the expansion is thus 
equal (up to a overall numerical factor) to 

d9 e^/^ exp[— TOT cosh — im — sinh6'l 



TVr — IX 



TVc + ix 



Ki {m\/ + x^Vc ^) 



exp [— TO V^T^ -I- X^Vc ^] 

\JVcT -t- ix 



(15) 



In the next step we use the fact that the single-electron 
Green's function factorizes into a charge and a spin piece. 
Let us consider intermediate states containing a single 
antiholon and any number of spinous. Factorization im- 
plies that if we carry out the sum over all multi-spinon 
contributions, we must obtain the conformal result ( pT| ) 
for the spin sector. Therefore we arrive at the following 
result for the two-point function 



(L,(x,r)Lt(o,0)) 



Zq exp[— m-\/ -|- x^Vc ^] 
2^V(«s 



T -h lX){VcT -h IX] 



(16) 



The leading corrections to (|T6|) involve intermediate 
states containing two antiholons and one holon and are 
thus of order 0(exp(— Smr)). Similarly we have 



(i?,(a;,r)i?t(0,0)) 



Zq exp[— to-\/t^ 



+ X'^Vc 



27r 



\J(vsT - ix){vcT - ix) 



(17) 



j- p^ were first written down by Wiegmann (for 
This remarkable result appears to have been 
subsequently forgotten and rediscovered only much later 
by Voitj— who conjectured the same form the Green's 
functionta. The. correct form was again reproduced by 
Starykh et alJtB, although on the basis of arguments 
which cannot be accepted as rigorous. An earlier attempt 
to calculate the singk-particle Green's function using the 
formfactor approachc2l led to an incorrect result yielding 
iCi-function instead of Ki in Eq.(p^). Finally, Parola 

and Sorella0 determined the spectral function by map- 
ping the problem of a single hole in a Mott insulator 
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onto an effective spin cfiain with skew boundary condi- 
tions. Their results are somewhat implicit so that it is 
difficult to compare them to (p7|). 

To conclude this section we would like to remark that 
the above results can be generalized to the case of the 
SU(N) Thirring model 



(L.(x,r)4(0,0)>^^^ 



1 



Stt {vsT + ix)i' 



VcT + ix 



Ki_±{m\/ + x'^i 



(18) 



Here we have used that the kinks in the SU(N) Thirring 
model have spin ^(l — 1-^)E3. A similar formula with 
X — > ~x holds for the right-moving fermions. 



A. Short-distance behaviour: RG calculations. 

Strictly speaking, single kink emissions in the charge 
sector allow one to consider only the frequencies up to 3m 
which is the threshold for the emission of two antiholons 
and one holon. However, since the numerical value of Zq 
is quite close to one eqns (|l^,|l^) are likely to provide 
a very good description of the Green's function in the 
entire frequency interval. Indeed, considering the short- 
distance (high energy) behaviour (m^/ + x'^Vc'^ <C 1) 
in (|l^) and (|l6|) we recover the standard expressions 
for the Green's functions of the Luttinger liquid (mod- 
ulo the prefactor Zq, which, as we have said, is al- 
most 1). A better description of the short-distance be- 
haviour is obtained by carrying out a renormalisation 
group (RG) analysis. The RG equations for the coupling 
constant g and wave function renormalisation Z in the 
case Vc — Vs = V are 



2 ^2 

50 -3+ |- In— +0(5') , 



Z = 1- 



(19) 
(20) 



where A is a UV cutoff and fi a subtraction point. The 
RG improved result for the two-point function is thus 



(i?,(x,T)i?t(0,0)) - 



Z{r) 



2t:{vt — ix) ' 



(21) 



where 



Z(r) ^ 1 - —g{r) , 
Ibn 



In(rm) ' 



r"^ = x"^ + v^T^ < m ^. 



(22) 



Comparing this to (17) for short distances, we find very 
good agreement. In the anisotropic case, where the short- 
distance behaviour of the Green's function is controlled 
by the Luttinger liquid exponent, the agreement becomes 



worse. For "large" anisotropy it may become necessary 
to take intermediate states with e.g. two solitons and one 
antisoliton into account in order to get a good description 
of the Green's function over a large range of energies. 



B. Retarded dynamical Green's function 

Given the expression (^) for the Green's function in 
position space we can calculate the dynamical Green's 
function by Fourier transformation and subsequent ana- 
lytic continuation to real frequencies. A straightforward 
calculation gives the following result for the retarded 
Green's function 



G(fi(^,g) = -Zo 



- + \/ m? + v^q^ — up' 



1 ~ a 
1 + a 



where a = Vg/vc- At a 1 this simplifies to 



UJ — vq 



m 



\J + v^q^ — up 



(23) 



(24) 



C. Spectral function 



The spectral function is defined as 



1. 



ARR(uj,q) = -;^ImG^(cj,g) . 

In the limit Vc — — v obtain from (| 
lowing simple result 



^i?J?,(t^,'7) 



Zq m Q{\uj\ — ^JiiP^^iPq^) 



Tx\u:-vq\ ^JZP 



(25) 
the fol- 

■ (26) 



An explicit expression for the general case Vg ^ Vc is 
easily obtained from (p^). We note that the result is 
exact for energies w < 3m as the intermediate states 
involving more than one antiholon will start contributing 
only at this energy. However, from previous experiencetj 
as well as on the basis of the arguments given above, we 
expect (23) to be an excellent approximation up to very 
high energies. 

For general Vc,Vs the spectral function is nonvanish- 
ing in a region with boundary uj = e{p) defined by the 
equation 



e{p) = inmq[vs\p - q\ + \J vlq^ + w?\ 
The solution of this equation is 



(27) 
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where 



VsP + mVl — 



Q 



eo(p) 



if P<Q 



ei{p) ifp>Q ' 



v/1 - i^s/v^y 



(28) 



(29) 



The spectral function generally has several kinds of sin- 
gularities. Firstly, there are singularities just above the 
threshold. For small momenta |g| ^ Q there is a square 
root singularity above the threshold uj = e(q) 



A{uj,q) 



Zq eo(g) + Vcq 
1 

X = 

y/uj- eo{q) 



(['7 + eo('7)][aeo(g) - q])' 



(30) 



There also is a square root singularity above the threshold 
for "large" momenta \q\ S> Q, where 



A(uj,q) 



Zo 



1 



1 



TT^/Vc - Vs y/\q\ -Q y/uj- ei(q) 



(31) 



This square root singularity is a general feature of the 
half-filled Mott insulator and apparently holds for any 
ratio of the gap to the bandwidth. Ga-lculations done 
for a model where this ratio is infinitec3 gave the same 
singularity. 

In addition there is a second square-root singularity for 
UJ — > eo{q) from below 



A{iO,q) 



Zq £o(g) + Vcq 
7T ^2eo{q) 
1 



{[q + eo(g)][wc9 - aeo{q)]y 



V^o{q) 



(32) 



For q ~ Q these two singularities collapse onto a single 
one and we have 



A{cj,Q) 



Zo fl 



1 — a 



2-U{Q)-{Lo-e{Q))-i (33) 



In Fig.|l| we show the spectral function for — Vg ~ 
u as a function of oj and q. As Q = cxd in this case 
the threshold is simply given by a; = m? -t- v'^q^ . The 
square root singularities above the threshold are visible 
as regions of high spectral weight. 




FIG. 1. Density plot of the spectral function j4_r_r(cj, g) as 
a function of uj and vq/m for Va = Vc = v. 

Fig.|^ shows a constant energy scan SlI oj — 3to. 

1 



^ 0.5 
S 




-4 -2 2 4 

vq/m 

FIG. 2. Constant energy (a; = 3m) scan of the spectral 
function for ct = 1. 

When we allow the spin and charge velocities to be dif- 
ferent an interesting effect occurs as can be seen from the 
density plot of Ajiji{uj,q) in Fig.^. For momenta q < Q 
the threshold occurs at the gap for a single soliton moving 
with momentum q and the spectral function looks very 
much like it did in the case a = 1. However, at momenta 
q > Q the threshold gets shifted to lu — Vsq + mVl — 
and a lot of spectral weight is concentrated between this 
threshold and the line ui — yjm?- + w^g^, where a second 
singularity occurs (see above) . This is shown in the con- 
stant energy scan Fig.^. For large momenta Vs,cq S> ra 
the double singularity at cy p-' Wc,s9 is similar to what 
occurs in a Luttinger liquidcj. 
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FIG. 3. Density plot of the spectral function Arr{lj, q) as 
a function of ui and Vcq/m for a — 0.4. 



9Z 

p{oj) ^ -2Im G(^) {uj, x = Q) = —6{\uj\ - m) . (34) 

VsVc 

The density of states (^) is constant for energies higher 
than the single-particle gap m. Formula (|4|) is exact for 
< 3m and expected to be an excellent approximation 
up to very high energies. 



IV. COHERENT SINGLE-PARTICLE MODE 
GENERATED BY THE INTERCHAIN 
TUNNELLING 

Taking the interchain tunneling into account in the 
RPA yields an expression for the single-electron function 
in the form of (0). The interchain tunneling gives rise 
to a branch of coherent excitations below the threshold 
of the ID spectral function A[ijj, q) in the region where 
T±(k) < 0. Therefore, even an infinitesimal interchain 
coupling leads to the coherent particle motion in the 
transverse direction and there is no confinement in the 
sense of Andersonn3. 



II 

3^ 
< 




FIG. 4. Constant energy (tj = 2m) scan of the spectral 
function for a = 0.4. 

The states located in the "wedge" ei{q) < uj < eo^q) 
correspond to a situation where the single antiholon car- 
ries only a small part of the total momentum and the 
remainder is made up by exciting (many) spinous. This 
effect is reminiscent of Cherenkov radiation produced by 
a heavy particle moving through a medium with a speed 
greater than the speed of light in this medium. 



D. Tunneling density of states 

From the Matsubara Green's function at coinciding 
space points in model (||) we extract the single-particle 
density of states 



A. Vs 



Let us first discuss the case Vs = Vc = v. For very small 
interchain tunnelling |Tj_(A;)| <C m the pole of (||) appears 
at a finite frequency very close to eoip) = \/i^'^-rn?. In 
the vicinity of the pole we have 



G{uj,q,k) 



Z{q,k) 
UJ - e(q,k) 



(35) 



where Z(q, k) varies strongly at small q and where for 
\vq\ « \m^/ZaTAk)\ 

/ \ 2 

ZQTi_{k)m \ 



e{q,k) = £0(9) - 



y/2eo{q)[eo{q) - vq] 



(36) 



In the chain direction, the dispersion of the coherent 
mode is asymmetric around q = and has a minimum at 
vq ~ [Zi:fT±{k)Y /m. For larger values of T±{k) this pic- 
ture remains qualitatively unchanged although the dis- 
persion law becomes more complicated. In the vicinity 
of the minimum the dispersion is approximately given by 



(g,fc) « Juj^ + v^{q- qoY 



(37) 



where wq, qo depend only on T±{k). At large positive 
vq ^ m we have 

eiq,k)^vq + ZoT^{k) , (38) 

and at large negative vq <ti — max(m, |ZoTi(A;)|) 



e(g, k) w -1 / + (vq)^ — [ 



ZoT±{k) -.2 
2vq -I 



(39) 
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FIG. 5. Dispersion of the coherent mode along the chain 
direction for several values of a; = ~ZoT±{k) and a = 1. 

The dispersion of the coherent mode depends on the 
transverse momentum only through Tj_{k). In Fig.^ we 
plot 



E,,iq) = e{q,k) 



(40) 



-ZoTj_(fc)=x 



as a function of q for several values of x. In Fig. 6 we 
plot the residue of the coherent mode for various x versus 
qv/m. 
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FIG. 6. Residue of the coherent mode as a function of 
the momentum along the chain direction for several values of 
X = -ZoT_L{k) and a = 1. 

We see that Z(q, k) is very small within the noninter- 
acting Fermi surface. 



B. Vs / Vc 

The case of different velocities can be dealt with in a 
completely analogous way. In Fig.^ we show the spec- 
tral function as a function of energy and momentum 
transfer along the chain direction for weak tunnelling 
ZoT±{k) = 0.5. The coherent mode is clearly visible 
below the threshold of the one dimensional spectral func- 
tion. 




FIG. 7. Density plot of the spectral function Arii{u!, q, k) 
as a function of lu and momentum transfer along the chain 
direction q for ZoT±{k) — 0.5. 

In Fig.|| we show a constant energy scan of the spectral 
function. 




vq/m 

FIG. 8. Aiir{uj, q, k) as a function of the momentum along 
the chain direction for u — 2m, a = 0.4 and ZoTj_(k) = 0.5. 

If we compare Fig.|| to the corresponding plot of the 
spectral function of uncoupled chains in Fig.^ we notice 
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that now there are two delta function peaks correspond- 
ing to the coherent mode, and that the singularities of 
the incoherent scattering continuum have been smoothed 
out. It turns out that most (~ 75%) of the spectral 
weight for fixed T±(k) (i.e. we integrate only over q) is 
located in the coherent mode at q2, about 23% sits in 
the incoherent spinon-antiholon continuum and only 2% 
is due to the coherent branch at qi . The sharp difference 
between the weights in the coherent modes is consistent 
with the picture presented on Fig. ^ 

As t± increases more and more of the spectral weight 
gets transferred to the coherent mode. For example, 
for ZoT±{k) = 3.2m approximately 95% of the spectral 
weight at energy w = 2m is located in the coherent mode 
at q2. However, it is important to note that the coher- 
ent modes dominate only a small portion of the Brillouin 
zone (the part where T± (k) is large) . If we consider the 
total spectral weight, i.e. integrate over the transverse 
momentum k as well, we find that the contribution of 
the coherent mode is generally small. 



V. FERMI SURFACE AND LUTTINGER'S 
THEOREM 

According to Luttinger's theorem the total number 
of particles in the system is proportional to the vol- 
ume of momentum space included in the surface defined 
by the singularities of lnG(w = 0,p). These singulari- 
ties may be due to either poles or zeroes of the single- 
particle Green's fpsction. As was recently emphasized 
by Dzyaloshinskiic^ the latter possibility must occur 
in Mott insulators thus removing an apparent violation 
of Luttinger's theorem by the interactions. Indeed the 
Green's function ( [2^ ) exhibits precisely this property: at 
Lu — it vanishes at q = 0, which corresponds to the 
Fermi surface of the non-interacting system. Since the 
position of the zeroes is unchanged in RPA, our results 
are in agreement with Luttinger's theorem. 

When t± exceeds a critical value, electron- and hole- 
like pockets of Fermi surface appear. The Fermi surface 
is determined by the equation 



Goi(0,g) = r^(k) 



(41) 



The volumes of electron and hole pockets are equal since 
rj_(k) = -T±{q + k) where Q = (7r,7r,0). So Lut- 
tinger's theorem continues to hold. 

Using the fact that at w = the Green's function is 
always real, we get from (p3|) 



^yvfq^ + Trl? — m 



sign(rj_(k)) = -sign(g). 



(42) 



The critical value of T± (k) necessary to produce the so- 
lution is 



V -fmi 



m 



= 3 



9a + I 



(43) 



where and is determined by the momentum q^ where 
the Fermi surface first appears 



3a - 1 + \/l + 10a + 9a2 
4a 



(44) 



The critical value of Zot""" varies from 2m at small Vg 
to ~ 3.3m for Vs ~ v^- 

The residue at the Fermi surface is given by 



2Zo jx-l) 
(1 + a) X 



{xa + 1) 



(x + 1) 



1/2 



(45) 



where x = -^/l + v'^q^ /m?. Near the critical value of 
Tl (k) the residue is numerically small but never goes to 
zero. For example, in the limit Us — > when a; = 2 we 
get 



Z{vs - 

and at Vs = Vc when x ; 

Z{vs = 



■ 0) « 0.58 Zo 
1.62 we have 
Vc) ~ 0.38 Zo. 



For a cubic lattice the Fermi surface forms an electron 
pocket around q- = go > 0, fcg = fcp = and hole pockets 
around q ~ —qo and fcg = ±7r, fcg — ±tt. The vol- 
ume of the electron pocket is the same as the sum of 
the volumes of the hole pockets. When the pockets are 
very small their shape can be determined by expanding 
Eq.(^2[) around the point {0,qo), using that 

T^ik)^T^iO)[l-kW] . 

Here kj^ denotes the deviation in the transverse direction 
from and 7 is proportional to the lattice spacing in the 
transverse directions a±. We obtain 



2TliO)j'kl + -iq~ qof - T^(0)2 - (tf^f , (46) 
where 



vl (xo + l)(4aa;o - 3a + 1) 



XQ{xa - 1) 



We also can estimate the anisotropy of the Fermi sur- 
face from (^). The anisotropy of the masses is given 

by 



m± 
mil 



vl (xo + l)(4aa;o - 3a + 1) 



xl{l + axo) 



(47) 



where we have further approximated T^{Q) « t™'". Us- 
ing that Vc — 2ia||, where a|| is the lattice spacing along 
the chains, we find that 
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^1- ^. A "I 

mil 7-^ 



(48) 



where ^ varies between 1.06 for Vg Vc and 0.62 for 
Vs —> 0. Thus the magnitude of the mass ratio is deter- 
mined by the competition of two factors one of which is 
large [t/m) and the other is small {a\\/^). As a result, 
the Fermi surface may not be very anisotropic. 

An obvious question is whether or not the RPA can 
be trusted to describe correctly the formation of a Fermi 
surface. An obvious problem of the RPA is that it au- 
tomatically reproduces a purely one dimensional result 
at the particular wave numbers p where the transverse 
hopping vanishes Tl(p) = 0. On a 2D square lattice this 
would be at the points py = ±7r/2. In the caae-pf coupled 
Luttinger liquids the improved treatment oilil indicates 
that in the vicinity of these points RPA is unreliable. In 
the case of coupled Mott insulators, the electron and hole 
pockets we find for sufficiently large t±_ are located at po- 
sitions far away from the points where 7±(k) vanishes. 

Recently a dynamical mean field theory approach has 
been developed, which replaces the quasi-lD system by a 
single effective chaixu Aom which electrons can hop to a 
self-consistent bathE3'Ea. The resulting model has to be 
analyzed by numerical methods and for coupled Hubbard 
chains it is found that for sufficiently large transverse 
hopping amplitudes an open Fermi surface (close to the 
one of the noninteracting model) is formed. This is in 
contrast to the electron and hole pockets we find in the 
RPA. 



VI. TRANSVERSE CONDUCTIVITY AND 
DENSITY OF STATES 

Using the results for uncoupled chains as well as the 
RPA expression for the Green's function of weakly cou- 
pled chains we can determine various other physical 
quantities. 



A. The transverse conductivity 

Let us consider a situation where the transverse hop- 
ping is only between nearest neighbour chains. The 
transverse current operator is then given by 

(49) 

where I is a chain index and x denotes the position along 
the chain direction. Using this expression we can an- 
alytically determine the leading contribution in t± to 
the transverse conductivity by using the result for the 
Green's function of uncoupled chains. We find 



2ZleH\ 1 

7r(wc - Vs) UJ 



arctan ■ 



(50) 



where 5 = {vc — Vs)/(vc + Vg). 

In the limit —i- ^ v this simplifies to 

= ^ -(2mMV(^/2m)2- 

TTV Ul 



1 



(51) 



We see that the transverse conductivity vanishes at 
the threshold u; = 2m and increases above it in a char- 
acteristic square root fashion. This is reminiscent ofLihe 
behaviour found for the longitudinal conductivity inEj. 



B. Density of States 

Within the RPA we can determine the density of states 
by integrating the RPA spectral function over all mo- 
menta. This needs to be done numerically. For simplicity 
we only consider the case^Up = Vg = v. For a 2D system 
with T±{ky) — tj_ cos(fcj,)c3 we obtain the results shown 
in Fig.|| As t±_ is increased, the Mott gap in the DOS is 
filled in and eventually a peak forms around zero energy. 



20 



15 



3, 



5 - 



0.2 



0.4 0.6 
CO/m 



0.8 



FIG. 9. Density of states for a 2D square lattice for fre- 
quencies smaller than the Mott gap. The curves are for t± — 2 
(solid), t± = 3 (dotted) and t± — 4 (dashed). 

The analogous analysis for a quasi 3D system with 
Ti_{ky, kz) = i± [cos(fcj,) -f- cos{kz)\ yields a qualitatively 
different result. As t± is increased the Mott gap is again 
filled in, but now the DOS around lj = Q remains basi- 
cally fiat and no peak develops in the regime where RPA 
applies, i.e. t±/m = 0{\). In order to understand this 
result it is instructive to consider the case 



7^^(k) 



D 

t± cos kj 



01. 



(52) 



In this limit ^^(k) can be considered to be a random 
variable with probability distribution 



Pit) 



1 



exp 



Dtl 



(53) 
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As we are interested in the formation of a Fermi surface 
we have to consider Tj_(k) = 0{m), which impHes that 
t±D/m — 0{1). This leads to the restriction that for 
large D 

« - , (54) 

so that the probabihty distribution (^ ) becomes ex- 
tremely narrow and only regions in k-space with Tj_ (k) « 
contribute to the DOS. However, in these regions there 
are no states at w ~ and there will therefore be no peak 
in the DOS at low energies. 

VII. BECHGAARD SALTS AS A POSSIBLE 
APPLICATION 

Our theory may be relevant to the Bechgaard salts 
andia particular to (TMTTF)2PF6 and (TMTSF)2PF6. 
(seec2l for a review). However, for our theory to be rel- 
evant, certain conditions must be met and this requires 
a discussion. The materials in question are quasi-one- 
dimensional and have a 3/4-filled band. The ratio of the 
hopping integrals in the three principal crystallographic 
directions is ta : ti, : tc = 1 : O-l : 0.005. Therefore, at 
sufficiently high temperatures one may neglect the hop- 
ping in the c-direction. Then one is left with a two- 
dimensional system of chains where each chain has only 
two nearest neighbours (A/" = 2). Such small number of 
nearest neighbours may put into question the applicabil- 
ity of RPA. This is, however, not the main worry. 

The principal problem one encounters in dealing with 
the Bechgaard salts is the problem identifying the correct 
low-energy effective theory. Let us neglect the interchain 
hopping for the time being and consider uncoupled ID 
chains. In the Bechgaard salts there are two separate 
mechanisms that have the potential to open a spectral 
gap. Firstly there are "double" Umklpjp processes due to 
the commensurate band filling 3/4oCj. These processes 
involve the scattering of four electrons at kp and four 
holes a.t —kp and are generated in the low-energy effec- 
tive theory by integrating out high energy modes. These 
processes will generate a gap only for strong interactions 
{Kc < 0.25). Secondly there is a small dimerization, 
which halves the Bdlkujin zone and gives rise to "single" 
Umklapp processesc3o. These open a gap already for 
weak interactions Kc < 1 but their coupling constant is 
proportional to the dimerization and thus small. At low 
energies the system is thus described by two indepen- 
dent Gaussian models - one for the charge and the other 
for the spin sector. The charge sector Gaussian model 
is perturbed by two operators: the 4fcF-harmonics (with 
kp being equal to 37r/2a) of the dimerization 

D ^ StY,i-irc++,^.Cn,a (55) 

n 

where St is a staggered component of the hopping in- 
tegral, and the 8kp component of the electron density 



operator (double-Umklapp processes). Since these oper- 
ators have different symmetry under parity transforma- 
tion (one is defined on links and the other on sites), in 
the continuous limit they are given by sin and cos respec- 
tively, such that the related Hamiltonian density is 

V = gi sin(v/87r7^c$c) + 92 cos(2v/8^$c) • (56) 

The behaviour of the system and applicability of our 
theory depend crucially on the value of K^. We can con- 
sider the following possibilities. |— in 

(i) The interactions are weak and Kc ~ 10c3. That 
the interactions are moderate is suggested by the renor- 
malization of the uniform magnetic susceptibility with 
respect to the value extraeted from the band structure 
calculations: Xs/xo ~ 2 — 3^. In this case, the (72-term in 
( p6| ) is irrelevant and the first term (due to dimerization) 
gives the sine-Gordon model which is equivalent to the 
charge sector of the Thirring model we have discussed. 
All our calculations are valid in this case. 

If we adopt this scenario we have problems with An- 
gle Resolved Photoemission (AREESj) data which do not 
show any traces of quasiparticlesO'Lj. We will return to 
the question about ARPES data later. 

(ii) There is a school of thought which advocates 
the small value Kc ~ 0.22 when both operators are 
relevanlEj. According to this school, the dimerization 
coupling is small and the g2-term dominates. This would 
agree qualitatively with ARPES data. 

If this point of view is correct, the calculations pre- 
sented in this paper are not applicable because, as it 
follows fronxJ, for such values of Kc the minimal form- 
factor corresponds to the emission of not one, but two 
solitons. We have discussed this scenario in our other 
publicationES. Here we cite just one conclusion from this 
paper: if Kc is small and the dimerization term is negli- 
gible, the value of the gap measured by ARPES should 
be equal to the optical gap observed in the frequency de- 
pendence of the on-chain conductivity, or twice the value 
of the activation gap in the temperature dependence of 
the dc conductivity: 

Aarpes = 2£:act = Aopt = 2to (57) 

Indeed, in (TMTTF)2PF6 ARPES measures a gap A « 
100 meV and thermal conductivitv-aneasurements give 
the activation gap i?act = 44meVEa, which is roughly 
one half. 

(iii) The third possible scenario is that Kc is small, but 
the dimerization is not negligible. To resolve this issue 
theoretically one has to estimate the bare values of the 
coupling constants 51,2- We know how to do this only for 
small U, where perturbation theory gives the following 
estimate 

52 - Eo{U/Eof . 

If both terms are relevant and the sign of 32 is positive, 
as it follows from the perturbative calculation, then the 
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gi-term leads to the confinement of the soUtons in the 
(72-sine-Gordon model. This is a rather difficult case for 
the theory because the double sine-Gordon model is not 
integrable. 

Let us assume for the time being that ifc ~ 1 
and see how our theory would square with experiment. 
(TMTTF)2PF6 is a Mott insulator; the Fermi energy is 
estimated as 115 meV, the hopping integral in the b- 
direction is of order of 14 meV, the optical gap Aopt = 2rn 
is of order of 900K. Here the transverse tunneling is not 
large enough to overcome the Mott gap. In trying to 
get a detailed comparison with our calculations one has 
to take into account that the ratio m/ep ~ 1/3 is not 
very small here, which reduces the chances of obtaining 
a quantitative agreement with any field theoretical ap- 
proach. 

(TMTSF)2PF6 is metaUic; the Fermi energy is esti- 
mated as 220 meV, the hopping integral in the 6-direction 
is of order of 20 meV, the optical gap Aopt = 2m is of 
order of 250K. This gives Tj_(0)/m « 4 so that the crite- 
rion for having a small Fermi surface ZQT±{0)/m > 3.3 
is satisfied here. Optical measurements for this mate- 
rial show a metallic Drude peak with a tiny amount of 
spectral weight (3 percent), separated by a gap from a 
very strong continuum. The metallic character pi these 
compounds is due to the transverse hopping^J^. For 
frequencies not too close to the gap the observed form of 
the optical conductivity in the direction along the chains 
is well described by the sine-Gordon modcE^. This fact 
together with the observed smallness of the Drude weight 
indicate that (TMTSF)2PFe is a good candidate for ap- 
plication of the present theory. 

An additional argument in favour of small pock- 
ets of Fermi liquid is that not all physical proper- 
ties of (TMTSF)2PF6 demonstrate the same degree of 
anisotropy. For example, the measured anisotropy of the 
plasma rfcequency for the Drude peaks is only of a fac- 
tor of 2c3. On the other hand, for the ratio of hopping 
integrals predicted by the band theory one should ex- 
pect it to be of order of (m^/m||)^/^ « 10. This fact 
in combination with the smallness of the Drude weight 
indicates that the Fermi surface is small and not very 
anisotropic. On the other hand, many of the properties 
of these material (especiaUjz. the magnetic ones) are typi- 
cally one-dimensional (seeEiJ for review) . Thus the overall 
picture is in reasonable agreement with the scenario we 
present. P-. 

The analysis of the Drude peak given inO indicates 
that the best fit can be obtained if one assumes frequency 
dependent effective mass and the scattering rates in the 
Drude formula: 



FiM = ro 



1 -t- a'^Lu-^ 



(58) 



cr(w) 



1 



47rri(w) 
m*{uj) ^ 

Wband 



• rm* (uj) 1 ' 
1UJ\ 

Aq 

l-ha2w2' 



This fit is rather suggestive because the frequency de- 
pendence is quadratic, like in Fermi liquid theory. This 
feature supports the point of view that the Drude peak 
comes from small pockets of Fermi liquid. 

There are other quasi-one-dimensional systems for 
which ARPES measurements have beerij)erformed: the 
blue bronze K0.3M0O3, which is metallic^, and the Mott 
insulator Sr2Cu03. The latter material, however, is not 
a good testing ground for our theory since the ratio of the 
Mott-Hubbard gap to the bandwidth is too large. The 
largeness of the gap precludes a detailed comparison with 
the results obtained in this paper. On the qualitative 
side the measurements demonstrate the appearencej-of 
two distinct dispersing maxima in the spectral functionEj, 
which is interpreted as a sign of spin-charge separation. 

ARPES is not the only way to ascertain the existence of 
quasi-particles. de Haas-van Alphen and Schubnikov-de 
Haas effects are perfect tools when one deals with closed 
Fermi surfaces. Obviously, the measurements should be 
made above the ordering temperature which may impose 
serious difficulties. 
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APPENDIX A: ELEMENTS OF RPA 

In this appendix we discuss some relevant aspects of 
the RPA in the interchain coupling. This is most easily 
done in position space. We denote the right and left- 
moving fcrmion operators by black and white circles. ID 
correlation functions are denoted by encircling a number 
of circles, the corresponding fermion operators are then 
all located on the same chain. In Fig.|l^ we show the 
corresponding diagrams for the ID two-point functions 
of right movers and left movers as well as the diagram 
fro the 2n-point function of right movers. Finally, we 
denote the interchain hopping matrix element ty (x — y) 
between sites x/a^ on chain i and site y/ao on chain j ^ i 
by a dashed line. We note that the hopping is local in 
time. 
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i,x i,y 



i,x i,y 



t:j(x-y) 



i,x j,y 



Gr...r(Xi,...,X2„) 



FIG. 10. Elements of the diagram technique in position 
space. Only one kind of many-particle "vertex" is shown. 

The first few diagrams in the expansion (in the inter- 
chain hopping) of the two-point function of right moving 
electrons for initial and final point located on the same 
chain is shown in Fig.|ll|. The contribution of a given 
diagram is obtained by summing over the positions of all 
"internal" circles, that is circles connected by a hopping 
line. 



(• • ) = ct 




+ 



FIG. 11. Diagrams for initial and final points on the same 
chain up to second order in the interchain hopping. 

The RPA expression for the (Fourier transform) of the 
single-particle Green's function is obtained by summing 
all diagrams that can be split into two parts by cutting 
any one hopping line, i.e. all diagrams of the type shown 
in Fig.pl 



FIG. 12. Diagrams extering the RPA expression for the 
single-particle Green's function. The sum is over all diagrams 
of the type shown. 

All diagrams neglected in RPA all contain loops. This 
enables us to embed RPA into a systematic perturbative 
expansion in a small parameter kq as follows. Let us 
consider an interchain hopping ^^(k) (for simplicity we 
take it independent of the wave number q along the chain 
direction) of the form shown in Fig.|l^ i.e. particle-hole 
symmetric but long ranged such that its Fourier trans- 
form is strongly peaked around the origin and the point 



(0,7r,7r). This means that in position space the hopping 
is long-ranged, i.e. the hopping amplitudes are of the 
same order within a range proportional to Kq^. 



FIG. 13. Schematic dependence of T±(k) on the trans- 
verse momenta q±. 



Every loop gives a contribution 
l5][T^(k)]V«^ 



(Al) 



and is thus suppressed. In this way one obtains a formal 
expansion in powers of Kq, the leading (k") order of which 
is given by the RPA. 



APPENDIX B: FORMFACTOR APPROACH IN 
INTEGRABLE QUANTUM FIELD THEORIES 

The spectrum of low-lying excitations of the half-filled 
Hubbard model consists of scattering states of gapped, 
spinless charge ±e excitations called holons and anti- 
holons and gapless, charge-neutcal excitations carrying 



spin ±i, the so-called spinonseS. We introduce labels 
h,h,s,s to distinguish between these four types of elemen- 
tary excitations. Their dispersions and exact scattering 
matrices areJcnown on the latticeC3 as well as in the field 
theory limitcS. As usual for particles with relativistic dis- 
persion it is useful to introduce rapidity variables 0c, s to 
parametrize energy and momentum 



Ea(Pc) — m cosh 0c ,Pa(Pc) — m sinli 
m 



P 



m 



-7 



a = h,h, 
= s,s. (Bl) 



Here we have set spin and charge velocities to 1. Let 
us now turn to the constuction of a basis of scatter- 
ing states of holons, antiholons and spinous. A conve- 
nient formalism to this end is is obtained in terms of the 
Zamolodchikov-Faddeev (ZF) algebra. The ZF algebra 
can be considered to be the extension of the algebra of 
creation and annihilation operators for free fermion or 
bosons to the case or interacting particles with factoriz- 
able scattering. The ZF algebra is based on thCpknowl- 
edge of the exact spectrum and scattering matrixcEl. For 
the SGM the ZF operators (and their hermitian conju- 
gates) satisfy the following algebra 
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zli9,)zl{9,) = zl,i92)zl,{e,)sti::koi - O2), 
z^^{e,)zl{e2) - z\{e2)siX^e, e,)z<{e,) 

+ {'^n)5ll5{ei~e2). (B2) 

Here S'^}''^H6) are the (factorizable) two-particle scat- 

tering matrices and Sj — s, s, h, h. 

Using the ZF generators a Fock space of states can be 
constructed as follows. The vacuum is defined by 



z,M\o) = 0- 



(B3) 



Multiparticle states are then obtained by acting with 
strings of creation operators Z^{9) on the vacuum 



|0„...ei)e„...ei =^t(e„)... 4(^1)10). 



(B4) 



In term of this basis the resolution of the identity is given 

by 



EE 

71=0 et 



dOi... dOn 

(27r)"n! 



\0. 



(B5) 



Inserting (B5) between operators in a 2-point function 
we obtain the following spectral representation 

71 



X |(O|O(O,O)|0„...0i), 



|2 



where 



f^{Oi . . . 9n)ei...e„ 

arc the form factors. 



(0|O(0,0)|^ 



(B6) 
(B7) 
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